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Arguably the most important problem in quantitative finance is to under- 
stand the nature of stochastic processes that underlie market dynamics. One 
aspect of the solution to this problem involves determining characteristics of 
the distribution of fiuctuations in returns. Empirical studies conducted over the 
last decade have reported that they are non-Gaussian, scale in time, and have 
power-law (or fat) tails However, because they use sliding interval 

methods of analysis, these studies implicitly assume that the underlying process 
has stationary increments. We explicitly show that this assumption is not valid 
for the Euro-Dollar exchange rate between 1999-2004. In addition, we find that 
fiuctuations in returns of the exchange rate are uncorrelated and scale as power- 
laws for certain time intervals during each day. This behavior is consistent with 
a diffusive process with a diffusion coefficient that depends both on the time and 
the price change. Within scaling regions, we find that sliding interval methods 
can generate fat-tailed distributions as an artifact, and that the type of scaling 
reported in many previous studies does not exist. 

Our analysis is conducted on one-minute intra-day prices of the Euro-Dollar exchange 
rate (obtained from Olsen and Associates, Ziirich) which is traded 24-hours a day. Let 
P{t) represent the exchange rate at time t and define the return of the exchange rate as 
x{t; t) = log [P{t + t)/ P{t)]. Here t represents a time during the day and r a time increment 
that is initiated at t. The analysis presented below is predicated on the assumption, for which 
we provide evidence, that the stochastic dynamics of x(r; t) is the same between trading days. 
Then, we find that the average movement taken over the approximately 1500 trading days 
during 1999-2004, (a;(r;t)) nearly vanishes for each value of t. A value of r = 10 min is 
used so that the autocorrelations in the signal P{t) have decayed sufficiently. The rest of 
our analysis is conducted on ffuctuations x(r; t) = x{t; t) — (x(r; t)) about the mean. 

A stochastic process has stationary increments if the distribution of x(r; t) is independent 
of t; otherwise, increments are nonstationary. Figure 1(a) shows the behavior of the standard 
deviation cr(r; t) = \/ (x(r; t)^) of the Euro-Dollar rate as a function of the time of day. If 
the stochastic increments are stationary, the curve would be ffat. Clearly, it is not. Instead 
(j(r; t) exhibits complicated nonstationary behavior while changing by more than a factor 
of 3 during the day. 

Our assumption of daily repetition of the stochastic process is validated by conducting a 
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corresponding analysis of fluctuations throughout a trading week p. Figure 1(b) shows the 
standard deviation of returns averaged over the 300 weeks studied. The approximate daily 
periodicity of cr(r;t) is evident, thereby justifying our approach. Similar observations were 
made on price increaments for Euro-Dollar rate in Ref. 6|. 

The standard deviation scales as power-laws with time during several intervals within 
the day. Power-law fits to the data in some of these intervals are shown by colored lines in 
Fig. 1(a). We focus our analysis on the time interval I which begins at 9:00 AM New York 
time and lasts approximately 3 hours. The data shown in red in Fig. 2(a) shows that the 
standard deviation within this interval scales like t'"^ where t is measured from the beginning 
of the interval and the index rj = 0.13 ±0.04. This scaling extends for more than 1.5 decades 
in time. Note that the value of rj is different for the other time intervals during which the 
standard deviation scales in time. Similar variation in scaling exponents during the day has 
been reported previously 0]. 

The scaling index within I does not change significantly during the six years studied. This 
is demonstrated by independently analyzing three two-year periods 1999-2000, 2001-2002, 
and 2003-2004. Figure 2(b) shows that the scaling index remains nearly unchanged between 
these two-year periods. 

We have also analyzed the behavior of other moments (|x(r; t)!'^)^/^ of the returns. Figure 
2(a) shows that each of the moments (3 = 0.5, 1.0, 2.0, and 3.0 also scales as a power-law 
in time, and furthermore that the scaling index for each of them is consistent with the value 
of 7] = 0.15. This nearly uniform scaling of the different moments suggests that the return 
distribution itself scales in time. Denote the distribution of x(r; t) by W{x,T]t), where 
the final argument reiterates that the distribution can depend on the starting time of the 
interval. In particular, when the increments are nonstationary W{x, r; t) depends on t. Our 
scaling anzatz is 

W{x,T;0) = ^nn) (1) 

where H is the scaling index, u = x/t^ the scaling variable and JF the scaling function. 
Note that the scaling anzatz is for a time interval starting from the beginning of I. 

In addition to scaling, the stochastic dynamics appears to have no memory. This can be 
demonstrated by evaluating the auto-correlation function 

. /, . , (x(t; ti)x{T;t2)) 
air: ti a r; to) 



We find that for r = 10, Ar{ti, ^2) = 1 if ti = ^2, and of the order of 10~^ when |ti — 12| > 10. 
This observation ehminates fractional Brownian motion [8] as a description for the underlying 
stochastic dynamics, and strongly indicates that dW{x,T;0)/dT depends only on x(r; 0) 
and T. If, in addition, W{x,t;0) has finite variance (see Fig. 4), it has been analytically 
established that the evolution of W{x, r; 0) is given by a diffusion equation 0, [l^ 

dW(x,T;0) 1 92 

g; =2 9^ iD{x, t)W{x, r; 0)) , (2) 

where D{x,t) is the diffusion coefficient. There is no drift term in Eq. (2) because x{T]t) 
has zero mean for all t. Note that the stochastic dynamics is completely determined by the 
diffusion coefficient, which, as shown below, depends on H. Hence, H can be considered to 
be the dynamical scaling index. 

Because we have found scaling, consider solutions of the form (1) to Eq. (2). When 
H = 1/2, the diffusion coefficient has been shown to be a function of u; i.e., D{x, r) = 'D{u) 
0]. If, in addition, ^'(m) is symmetric in -u, it is related to the scaling function by T{u) = 
D(u)~^ exp (— /" dy v/D(v)) 0, [ill. When H 7^ 1/2, we can "rescale" time intervals by 
f = In f, the stochastic process has a scaling index 1/2 and a diffusion coefficient 

of the form V{x/^/f). Converting back to r, D{x, r) = 2Ht'^"-^V{u) Q. 

Statistical analyses of financial markets have often been conducted using sliding interval 
.e..od. BflQfiflQQ. a..n,e t.a. a. .a«on„. even 

if they are not. For example, they compute the distribution Ws{x, r) = {W{x, r; t))t, where 
{.)t indicates an average over t. Many of these studies have reported that Ws{x, r) scales as 

Wsix,T) = -^Mv), (3) 
where v = x/t^^ and 1/2. It has also been reported that the scaling function JF^ has 



power-law (or fat) tails |4, |^. However, it is important to understand that Ws{x,t) is a 
solution of Eq. (2) only when the stochastic process has stationary increments, in which case 
H = Hs = 1/2. In general, Hs and Ws{x, r) are different from H and W{x, r; 0). Next, we 
give an explicit example where this is the case, and, in addition, Ws{x,t) appears to have 
fat-tails even though W{x,t; 0) does not. 

Consider a diffusive process initiated at x = that has a variable diffusion coefficient 
2Ht'^^^^{1 + \ u\). Its distribution has a scaling index H and a scaling function J^{u) = 
^exp(— 0, 11|. (See the discussion following Eq. (2).) Numerical integration of the 



stochastic process for H = 0.35 confirms this claim, see Fig. 3(a). In contrast, Ws{x,t) 
calculated from the same data appears to scale with an index Hs = |. Unlike JF which is bi- 
exponential, the apparent scaling function J-'s (shown in Fig. 3(b)) has fat-tails. However, a 
careful analysis reveals that distributions Ws{x, r) do not scale in the tail region, and hence 
that J-'s is not well-defined. Differences analogous to those between H and Hs have been 
noted for Levy processes |3] and for the R/S analysis of Tsallis distributions isl ]. 

The behavior of cr(r; t) (Fig. 2(a)) can be calculated for variable diffusion processes. 
Assuming that r is small, Ito calculus gives = x(r;t)^ = D{x,t)T. Averaging over 
returns at t gives 



(4) 



dxW{x,t; 0)D{x,t) 

In a variable diffusion process, W{x,t;0) = t~^J^{u) and D{x;t) = 2Ht^^~^V{u)] conse- 
quently 



(5) 

independent of the exact form of T>{u). Results for the Euro-Dollar rate within the interval 
I (Fig. 2(a)) which showed that r] ^ 0.15 are therefore consistent with a scaling index 
H = ^ — 1] ^ 0.35. Note that, unlike for Levy processes and fractional Brownian motion, 
H < 1/2, and is substantially less than Hs reported in previous analyses of the Euro-Dollar 
exchange rate (between 0.5 and 0.6) ^,14, 31. A general calculation for the moments of a 
variable diffusion process gives 



{\6xfy/^ r^t""-'/', (6) 

for all (3, consistent with results shown in Fig. 2(a). 

In order to estimate Hs for an arbitrary variable diffusion process, we note first that 
{x{t + r; 0)^) = {x{t; 0)^) -|- (a;(r; t)^) for any diffusive process without memory (see Ref.0|). 
Then, using the scaling anzatz (1), setting c = j du uJ^{u), and taking the sliding interval 
average 

(x(r; tf), = (c(t + r)2^ - ^ 2Hc{t'''-')tT, (7) 

where the last approximation is valid when r ^ t, a condition that is true for most intervals 
of length r in a sliding interval calculation. Hence (x(r;t)^)( ~ r. Consequently, Hs = 1/2 
regardless of the value of H\ 
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Finally, we introduce a method to extract the empirical scaling function T from the 
Euro-Dollar time series. Unfortunately, the available data are insufficient to determine T{u) 
accurately using the usual method of collapsing W[x^ r; 0) for multiple values of r. However, 
since we have determined 0.35) independently, we can use Eq. ^ for multiple values of 
r in the interval I (i.e., r between approximately 10 and 160 minutes) to determine T . The 
result is shown in Fig. 4(a). Note that the distribution has an approximate bi-exponential 
form. Since exponential distributions have finite variance, all assumptions needed for the 
derivation of Eq. Q are justified. However, it is asymmetric and decays more slowly on 
the negative side. By contrast, the empirical sliding interval scaling function Tsip) for the 
same time interval is shown in Fig. 4(b). For this case, the scaling collapse is achieved for 

nn 

Hs = 1/2. Fs{v) appears to have fat tails, consistent with previous reports |5|,ll6[. However, 
in light of the example discussed earlier and the fact that 7^ 1/2, it is unlikely that JF5 is 
well-defined for this financial market data within the interval I. 

Variable diffusion processes exhibit another signature (stylized fact) of market fluctua- 
tions. Although their autocorrelation vanishes, a large fluctuation will typically produce a 
large value of and hence a return with a large diffusion coefficient. Consequently, a large 
fiuctuation is likely to be followed by additional large fiuctuations whose signs are uncor- 
related to the first |,9j. As a result, the autocorrelation function for the signal |x(r;t)| (or 
for the signal x(r;t)^) will decay slowly in t. Such behavior, referred to as the "clustering 
of volatility" is seen in the Euro-Do 



studies of other financial markets 
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lar exchange rate and has been reported in empirical 

The analysis given here applies to stochastic dynamics of a single scaling interval. How- 
ever, the daily fiuctuations in the Euro-Dollar rate are a combination of scaling intervals 
with distinct scaling indices, and possibly regions with no scaling. We have not yet deter- 
mined how to extend our analysis beyond a single scaling region. Bacuase of this, it is not 
clear how to interpret the distributions over intervals longer than a scaling region, including 
inter-day data. 

We have shown that stochastic fiuctuations in the Euro-Dollar rate have uncorrelated 
nonstationary increments during the course of a trading day, and that there are intervals 
during which their absolute moments scale like a power-law in time. The stochastic dynamics 
during these scaling intervals can be described by a diffusion process with variable diffusion 
coefficient. We have also shown that sliding interval analysis of variable diffusion processes 
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can give an incorrect scaling exponent and in addition can give scaling functions with fat- 
tails even when the underlying dynamics do not have them. Indeed, this appears to be the 
case within the interval I. 

The authors would like to thank A. A. Alejandro-Quinones for discussions. They also 
acknowledge support from the Institute for Space Science Operations (KEB, GHG) and 
the NSF through grants DMR-0406323 (KEB), DMR-0427938 (KEB) and DMS-0607345 
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FIG. 1: (a) The standard deviation (T(r;t) = \J {x{t; t^) of the daily Euro-Dollar exchange as a 
function of the time of day (in GMT). Here r = 10 min to ensure that autocorrelations in P{t) 
have decayed sufficiently. Our statistical analysis assumes that x{T;t) follows the same stochastic 
process each trading day. The average indicated by the brackets (.) is taken over the approximately 
1500 trading days between 1999-2004, and the standard error at each point is typically 3%. Note 
that, if the stochastic dynamics had stationary increments, (t(t; t) would be constant. Instead, it 
varies by more than a factor of 3 during the day, thus showing explicitly that the exchange rate has 
nonstationary increments. Notice also that cj(r; t) scales in time during several intervals, four of 
which are highlighted by colored lines that are power-law fits. Our analysis focuses on the interval 
I shown by the horizontal solid line, (b) The weekly behavior of cj(r; t) for the same data. Observe 
that it exhibits an approximate daily periodicity, thereby justifying our assumption of the daily 
repeatability of the stochastic process underlying the Euro-Dollar exchange rate. 
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FIG. 2: (a) A log-log plot of {x{T;tf)'^/l^ for /5 = 0.5,1.0,2.0, and 3.0, demonstrating power law 
decay i"'' for each index. Here t is measured in local New York time stating at 9:00AM. The data 
for (3 = 0.5, 1.0, 2.0, and 3.0, shown in blue, green, red, and black, respectively, have scaling indices 
(given by the slopes of the solid lines) r/ = 0.15 ±0.02, 0.14 ±0.02, 0.13 ±0.04 and 0.13 ±0.08. All of 
these values are consistent with r] ~ 0.15, and hence a dynamical scaling index of H = 1 — 0.35. 
The error estimates on the exponents are the standard errors from the nonlinear fit including the 
standard deviations for each time point, but neglecting any correlations between them, (b) The 
behavior of the standard deviation cj(r; t) in the interval I during each of the periods 1999-2000 
(blue), 2001-2002 (red), and 2003-2004 (green). The scaling index from nonlinear fits for the three 
data sets are 0.13 ±0.06, 0.14 ±0.04 and 0.14 ±0.07. The near equality of these indices shows that 
the scaling index is nearly invariant over time. 
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FIG. 3: (a) The scaling function of the return distribution T calculated from a collapse of data 
for r = 10 (blue), 100 (green), and 1000 (red) units. The results are from a set of 5,000,000 
independent stochastic processes with variable diffusion. The scaling index used was H = 0.35, 
and the diffusion coefficient was 2Ht^^~^{l + \u\). Note that J- is bi-exponential, as discussed 
in the text, (b) The sliding interval scaling function J^s calculated from the same runs. Shown 
are results for sliding intervals with r = 10 (blue), 100 (green) and 1000 (red) units from runs of 
length 10, 000 units. Unlike J^, it appears to have fat tails. The scaling index used here for which 
the scaling collapse is achieved is Hs = 1/2 even though the dynamical scaling index is H = 0.35. 
Note, however, although the central part of the distribution scales well, the tails do not. 
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FIG. 4: (a) The empirical scaling function T for interval I calculated assuming the scaling anzatz 
Eq. (1) with H = 0.35 and values of r between 10 and 160 minutes. Note that T is slightly 
asymmetric and approximately bi-exponential. Since exponential distributions have finite variance, 
all assumptions needed for the derivation of Eq. @ are justified, (b) The empirical sliding interval 
scaling function Ts for interval I calculated by scaling collapse of data using the anzatz Eq. (3) 
for r of 10 (blue), 20 (green) and 40 (red) minutes. Note that Ts has fat-tails. 
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